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Abstract 

We focus on the definition of the unitary transformation leading to an effective second order 
Hamiltonian, inside degenerate eigensubspaces of the non-perturbed Hamiltonian. We shall prove, 
by working out in detail the Su-Schrieffer-Heeger Hamiltonian case, that the presence of degene- 
rate states, including fermions and bosons, which might seemingly pose an obstacle towards the 
determination of such " Frdhlich-transformed" Hamiltonian, in fact does not: we explicitly show 
how degenerate states may be harmlessly included in the treatment, as they contribute with va- 
nishing matrix elements to the effective Hamiltonian matrix. In such a way, one can use without 
difficulty the eigenvalues of the effective Hamiltonian to describe the renormalized energies of the 
real excitations in the interacting system. Our argument applies also to few-body systems where 
one may not invoke the thermodynamic limit to get rid of the "dangerous" perturbation terms. 
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I. INTRODUCTION. 



Effective Hamiltonians obtained by unitary transformations truncated to second order in 
the perturbation, which here we shall label as of "Frohlich-type" , are ubiquitous in condensed 
matter physics 1 . Given a starting Hamiltonian H — H + I, with a basic part H and a 
perturbation / , one may write the unitary transformed Hamiltonian neglecting terms of 
third or higher order in the perturbation as: 

H = e n (H + I)e- n = H + / + [TZ, H ] + [TZ, I] + \ [TZ, [TZ, H }} + O (i 3 ) (1) 

The antihermitian generator TZ of the unitary transformation e 11 is determined by im- 
posing the vanishing of the first-order terms 

I+[K,H o ] = (2) 

so that the transformed second-order Hamiltonian reads 

H {2) = H + ^[n,I] + O{l 3 ) (3) 

To solve Eq.2 for TZ , one may explicitate it on a complete set of orthonormal eigenstates 
\X a ) of Ho : 

H \X a ) = S a \X a ) l = J2\X a )(X a \ (4) 

a 

yielding the basic equation to determine the matrix elements of TZ 

(X p \I\X a )=(E p -E a )(X p \Tl\X a ) (5) 

In solving Eq.5 one has to distinguish between the cases of non-degenerate or degenerate 
eigenstates. For non-degenerate cases, the matrix elements of TZ are directly obtained and 
the second-order Hamiltonian is straightforwardly derived, as usually done. The presence 
of degenerate eigenstates of H requires a special analysis. Namely, if £ a = £g the right 
hand side of the equation vanishes identically for any finite {X^\TZ\X a ) , and two possibilities 
appear, depending on the perturbation /. If also (Xp\I\X a ) = , the left-hand side vanishes 
as well, and one gets no condition for the corresponding matrix element of TZ between these 
degenerate states. On the other hand, for (Xp\I\X a ) ^ 0, Eq.5 would make no sense except 
for the case of a divergent (Xp\TZ\X a ). We will focus on this problem, and explicitly show 
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by working out in detail a non-trivial example that a well defined unitary transformation 
can nonetheless be defined when degenerate eigenstates are present (as first observed in 
general by Wagner 2 ) , and through it the second-order Hamiltonian be obtained. To fix 
ideas, we will consider the case of a non-perturbed Hamiltonian H describing a system of 
non-interacting bosons and fermions, in particular the Su-Schrieffer-Heeger electron-phonon 
Hamiltonian. The degenerate states \X a ) and \Xp) might differ in having one more, or one 
less, boson, which is created or destroyed in a real (as opposed to virtual) scattering process 
with the fermions, as described by perturbation J, conserving the total energy between initial 
and final states. We show how the generator 1Z can be defined, such that the frequently 
used procedure of identifying renormalized excitation energies of the real interacting bosons 
with eigenvalues of their effective Hamiltonian, obtained by averaging over the fermion 
wavefunction, can be applied without difficulties. 



II. THE SU-SCHRIEFFER-HEEGER ELECTRON-PHONON HAMILTONIAN 



Let us work out in detail a specific example based on the Su-Schrieffer-Heeger electron- 
phonon Hamiltonian 1 , for which the non-interacting Hamiltonian reads: 

ka q ^ ' 

Introducing the number operators n ka = c\. a c ka , for itinerant electrons in states character- 
ized by: crystal momentum k , spin a, and energy E k , and v q = b q b q , for phonons in states 
characterized by crystal momentum q and energy Q q , the generic eigenstate \X a ) = \X k ) 
of H may be written as a product of fermonic and bosonic orthonormal eigenstates: 



Xkqo) — \-- n kai n 



k+q,ai n k-q,ai 



).\..V„, V _ 



q- 



Above, we have explicitly written only the occupation numbers of the states of interest for 
our discussion. For instance, the boson state indicated simply as \v q --), denotes the following 
N -phonon state in occupation number representation (N different wavevectors) : 



(bir\o) 



while the orthonormal boson states ({v q \v p ) = 5 pq ) are such that: 



(7) 



(8) 



The exact eigenenergy S kqa of \X kqa ) is: 

Zkqa = E k* ("far) + E k+q,a ("fc+?,<r) + E k-q,a {n k _ q ^) + Ml q (u q ) 

+ E EpvM+Yt^M (9) 

p=^(k,k±q) r^q 

For this Hamiltonian, the perturbation term / reads: 

1 = -jfr E r *,*- ? c U- ? ,. (&L S + 6 S ) (10) 



kqa 

where the bond-stretching interaction amplitude, resulting from a modulation of the electron 
hopping tij , for a lattice with a centre of inversion is : 

V k,k- q = * E 9ij ( sin [(* - 9) A y] - sin [ fcA y] } % = dift-Rj) 

Notice r fejfc _ ? = -r* !fe _ 9 . Thus, rewriting £ feqCT T^^c^c^^ = 
Efe,* T k,k~q c l- q ,a c k,* b i > we can decompose the perturbation as 

7 = -4E( 7 *V+ 7 ^) (12) 

v fegcr 

where 

^fcgcr = ^k,k-q C k-q,a C k,a^q ^kqa = ^ k,k-q C ka C k-q,a^q (13) 

A given perturbation term I kq a yields non- vanishing results only when applied to the 
specific state \A kqcr ) where: 

\ A kqa) = Iha^k+q^^k-q^-M^q^-q-) \ A kqa) = I ka, l k-q,*i •> I "-q, "q- •) ( 14 ) 

Namely, 

= T lk- q \/^ q ) + l|0 te , W>K + 1, *- ff > = n,k-q\K°) ( 15 ) 

4^1^) = r^-gV^I 1 ^'^-' ^-)'^ ~ l ' V ~q) = T k,k-q\ B kqa) ( 16 ) 

where we have introduced the notation \B kqa ) for the states resulting from applying the 
perturbation. Notice also that: I^g J A k q a) = 0> anc ^ that spates ) and \B kqa ) are 
orthogonal: 

(^feqal-^feqa) = fi ( A kqa\^kqa\ A kqa) = ( 17 ) 

To fix ideas, in the following we will solve in detail the problem for the case of two degenerate 
eigenstates of the unperturbed Hamiltonian. The argumentation can be straightforwardly 
extended for larger size of the degenerate eigensubspace, without altering the conclusions. 



III. THE CASE OF TWO DEGENERATE STATES. 



Let us now assume that the spectrum of non-interacting energies of the system is such 
that, given the state \A+ q(r ) = \l kr7 ,n k+qa ,O k _ q(7 )\v q ,i;_ q ) with energy E£ qtr , a phonon exists 
such that \A+ qa ) is degenerate with the state \B+ qa ) = y/(u q ) + l\O ka ,n k+qa , l k _ gt0 )\u q + 
1, v_ q ) . Thus £ kqa = £ kqa , which implies E k = E k _ q + Ml q . In this case, we might be in 
trouble with Eq.5, as mentioned in the Introduction. 

To tackle the problem, it is convenient to rewrite the condition for the transformation 
generator 1Z , Eq.2 , specifically isolating in I the "dangerous" terms I kqa + I kqa : 

1 = -7^ E Wr + (i - w.) + -4 K« + 4v) = J i + J ^ (is) 

v p rT V 

Analogously, we write for the generator: TZ = 1Z\ + lZ kqa , explicitating its restric- 
tion to the degenerate subspace in last term. Notice that the commutators [Ii,I kqa ], 
[IZi, IZkqa] , [R-i, Ikqa] and [Ii,TZ kqa ] vanish, since the operators to be commuted involve 
different wavevectors. Thus, the transformed Hamiltonian of Eq. 1 now reads: 

H = e ( ni+nk «°) (H + h + l kqa )e~( ni+nk «°) = 



= H + I x + l kqa + [{IZi + 1Z kqa ) , H ] + [1Z 1: h] + [TZ kqa ,l kqo .} 

+i [n u [H U H ]\ + ^ [JZ kqa ,H ]] + ^ [1Z kqa , [n u H ]] + ^ [Kkga, [JZ kqa ,H }}+O (/ 3 ) (19) 

The condition to be imposed in order to eliminate the terms linear in the perturbation 
takes the form: 

h + T kqa + [(fti + K kqa ) , H Q ] = (20) 
The above equation for 1Z can further be split into two independent constraints 

h + [fti, ff ] = [fti, if ] = -h (21) 

^fcgcr + [T^fcgtr, -f^o] — [R-kqa, H ] = —T kqa (22) 

Therefore, using these two constraints in Eq.19 and dropping all vanishing commutators, 
it follows that to second order one has: 

H^ = H + l - [K u h] + \ [K kqa , Z kqa \ + O (J 3 ) (23) 

where the matrix elements of commutator [7^i,ii] can be evaluated safely even over the 
pair of degenerate states. In the following, we shall focus on the matrix elements of the 
potentially "dangerous" commutator \Ti kqa ,X kq(J ] . 



A. The action of the perturbation on the degenerate states. 



As \B£ ) is not normalized, let's introduce the corresponding normalized eigenstate \B kqa ) 



as: 



Bkqa) ~ n v -, \Bkqa) ~ \®kai n k+q,ai ^k-q^Wq + 1; V -q) 

V\ V q) + 1 



(24) 



Let us evaluate the matrix elements of (l kqa + I kqa ) inside the degenerate sub- 
space spanned by: \A kqa ) : \B kqa ) .The diagonal elements vanish by orthogonality: 
( A t q J Ka + Kqa) \ A tqa) = (^tqJ {Uq* + Kqa) \®tq«) = since the perturbation changes 
the number of bosons, while: 



(Btqal (ikqa + hqo) \ A kqa) ~ ^k 



k,k-q\/ \ V ql 



+ ^=[( A tqj(lL + Ikq„)K«)Y (25) 



B. The perturbation-split equivalent states and the corresponding matrix ele- 
ments. 



By a rotation inside the degenerate subspace, we now diagonalize the restriction of the 
perturbation^ ( l£ qa + I^ qa ), there. From its eigenvalues, we obtain the respective eigenen- 
ergies of the Hamiltonian including first-order perturbative corrections: 



p± _ cA I 
^kqa ~ °kqa ^ 



i^kqal /Jy {Hqa + hqa) \ A kqu) 



£t ± ^\ V Kk~q\\](Vq) + ^ (26) 



with respective orthonormal eigenvectors of Hq : 

ivc> = ^ W4W + Ka)} \r kqa ) = ^ [MAtj - \B + kqa )] 



A EE 



^k,k-q 

|r 

I k,k— q\ 



i\' where A' = sgn 



Y,9u {sin [(k - q) A y ] - sin [A; Ay] } 



(27) 



(28) 



The diagonal matrix elements of "dangerous" perturbation terms in the \^ kqa ) subspace 



are: 



(rtqJ [!tqc + ^qa) \^tqa) = ~ | F \,k-q\ \j <^> + 1 
WkqJ (4V + J kqa) l^kqa) = \ T k,k~q\ \j W + 1 

while the off-diagonal ones vanish. 



(29) 
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C. Matrix elements of the effective Hamiltonian in the \ip^ qa ) , \ipk qa ) subspace. 

Of interest is the effect that the " dangerous" perturbation term might have on the matrix 
elements of of Eq.23. We will show that, with an appropriate definition of 1Z, all matrix 
elements of \lZ kqcr ,T kqa \ in the subspace spanned by \ip kqa ) and \ip kqa ) vanish. 

For the evaluation of the off-diagonal elements of [lZ kqa .,Tk qo ] it will be convenient to 
use the following decomposition of the identity operator: 

1 = E \ X pro)( X pra\ i 1 ~ ( A pra\ X pra)) (l - (B+ ra \X pra )) 
prcr 

+\i>t q *)(ri q *\ + \^ q *)(^ q *\ (30) 

where use has been made of the equality: I^X^I + I^V> (faqJ = \ A k q a)( A k Q a\ + 
\^kqa)(^tqJ- ^ OT simplicity, in the following we will assume that except for the subspace 
explicitated above, H possesses no other degenerate eigenstates. 
We have: 

i^kqal [^kqcr^kqcr] \^qa) = 

I 

= "y= [(^kqJ^kqa (4qa + ^qa) l^kqa) ~ (^qal ^kqa + ^qa) ^kqaWkqa)] ( 31 ) 

Inserting the identity operator 1, the first contribution to the rhs of Eq.31 , reads: 

i^tqjKkqa i 1 ) ^kqa + ^qcr) l^kqj = 

= E i 1 ~ ( A U\ X pra)) (1 " (K<r\ X P ra)) (^tqal^kqa \ X pra) ( X pra\ + \^kqa) + 

prcr 

+ ^tqa\' R kqa\^tqa)^tqa\ { J kqa + 4^) \^kqa) + k^kqc l^kqcr l^kqa) (^kqa I { J kqa + ^qcr) \^kqa) 

(32) 

For the first term, notice that J2 prC r involves only eigenstates outside the degenerate 
subspace specified, while, using Eqs.15 , 24 and 25, one has: 

( X pro\ {Hqa + hqa) l^kqa) = 
= -J=r*k,k-q\l (Vq) + 1 ( X pra\Bkqa) ~ k,k-q\J ( U q) + 1 ( X pra I A tqa) = ( 33 ) 

where the last equality follows from orthogonality, in both scalar products. 
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The next term in Eq.32 vanishes because {^> kqa \ (lk qa + hqo) I^Wr) = • The last term 
in Eq.32 , due to Eq. 29, reads: 

i^kqal^kqal^kqj^kqj {^qa + ^qa) l^kqa) = 
= ~ Kk- q \ 7^) + ! (^kqal^kqa^kqa) (34) 

The off-diagonal matrix elements of lZ k are required to satisfy the condition given 
by Eq.5. But due to the diagonalization of the perturbation performed, and H .\ip^ ) = 
^kqal^kqa) > one * s m the case m which of both sides of Eq.5 are identically zero. Thus, one 
in fact has the freedom to choose the value of {^tqaO^kqaWkqa) '■ as zero > m particular. In 

the following, we will justify this choice as a reasonable and most convenient one. 

First, one can justify that choice by an analytical continuity argument. Having performed 
the diagonalization of the perturbation, {4> kqa \ {l^qa + ^kqa) li'kqa) = ® holds exactly for any 
value of |l\ fc _J > 0. If in Eq.5 one would replace the zero-order (degenerate, here) 
eigenenergies difference by the vanishing perturbation limit of the first-order eigenenergies 
splitting, 

= |rH0 ( E * qa ~~ E ^ ^kqal^kqal^kqa) = 
= ir^O^ l r *'*-«l + 1 ^tqa^kqMkq.) (\* ' k,k-q\ > 0) (35) 

one would see that in the whole neighbourhood of |r fc /,_J = the consistent definition of 
the transformation would be: 

(ri qa \K kqa \^ qa ) = (\T ktk _ q \ > 0) (36) 

By continuity, it thus seems reasonable to choose (i'tqcX^-kqa^kqu) = 0> a ^ so f° r l^^-J = 

0. 

Further, as a hand- waving argument, one could also mention that since TZ k has the same 
operatorial structure as (^ (T + -^ (7 ), one might expect the two objects to have similarly 
vanishing off-diagonal matrix elements in the \ip kq(7 ) subspace. 

Going back to Eq.31, with our choice we thus have a vanishing contribution for the first 
term of the commutator. To evaluate the second term of the commutator we again introduce 
the identity decomposition previously used, yielding: 

(i>kq*\ (4V + tkqo) (!) ^kqal^kqa) = 

8 



pro- 

+ (^«xl (4 + go + 4^) l^<x)(^<xl^l^) + (^1 (Ikqa + Jkqa) \^kqa) i^kqal^kqal^kqa) 

(37) 

Here, the first term is zero by orthogonality, analogously to Eq.33, and having chosen to 
define: {fpkqal^kqal^kqa) = 0' the secon d term of Eq.37 vanishes as well. Finally, the last 
term also vanishes, because \^ kqcr ) are orthogonal eigenvectors of the perturbation. 

Therefore, Eq.31 yields vanishing off-diagonal matrix elements of \lZ kqa ,Ik qa \ ■ 

Let us now evaluate its diagonal elements. We have to consider 

= -1= l(i>ijn kqa (/+„ + i kq(7 ) \^ qa ) - (^j (/+, + ^ CT |V£ CT >] (38) 

Again, inserting the identity decomposition, the first contribution to the commutator in 
Eq.38 reads: 

i^qjKkqa ^kga + ^qa) \^tqa) = (^tqal^kqa (l) (^kqa + ^qa) \^tqa) = 

= (l ~~ (A pra \X pra )) (l — (B pra \X pra )^ (lp£ qa .\1Z kq(J \X pro .) {X pra \ (l kqa + Ikqa) IV ; fcgo-) + 

pru 

+ {^tq„\K kq Mt qlT ){^tqc\ (tkqa + ^qa) \*l>kqa) 
+ ^W R kq^kqJ^kqJ {^q* + 4^) 1^) = (39) 

The last equality follows since in the first term the matrix elements (X pra \ {l kqa + I kqa ) \iPkqo) 
vanish by orthogonality condition, as we have shown when evaluating the off-diagonal matrix 
elements; in the second term the diagonal element {^tqcX^-kqaWtqa) = ® > due to an tiher- 
miticity of TZ kqa ; and, in the last contribution, the perturbation diagonalization yields 

(iKi+u i^) = o. 

Analogously, it is immediate to demonstrate that the second term of the commutator 
also vanishes, leading to: 

+ Ikqa)] K.) = (40) 

Using the same arguments, one can easily show that also the other diagonal element 
vanishes: {^J [K kq(T , + I kqa )] |^> = 0. 
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In conclusion, all contributions of the " dangerous " perturbation terms to the second-order 
transformed Hamiltonian have vanishing matrix elements inside the degenerate subspace, 
if one avails oneself of the freedom one has to complete the definition of generator 1Z and 
chooses that its restriction to the degenerate subspace is diagonal in the same basis of 
eigenvectors of the perturbation, there. This conclusion is clearly not dependent on the size 
of the degenerate eigensubspace. In our example, the choice is: {^tqcX^kqaWkqa) = 0> wnere 
l^kqo) diagonalize the perturbation restriction to the two-dimensional degenerate subspace. 
Therefore, the eigenvalues of the complete Hamiltonian in the degenerate subspace only 
have contributions from the "innocuous" commutator [i?i,7i]. Using Eq.23 , we can indeed 
write: 

<CJ# (2) IVC> = <^VI {#o + 1 [Ri, h] + ^= [n kqa , (i+ a + r kq(J )} } < CT > = 

= (^Vl {^o + \[Ri, h] } |< CT > = S A 5 afS + [Ri, h] |< CT > (a, = ±) (41) 

IV. DIAGONALIZATION OF ll ' 1 IN THE DEGENERATE \^ qa ) SUBSPACE. 

At this stage, to use the "Frohlich" -transformed Hamiltonian in the presence of degenerate 
states of Hq it would seem necessary to actually determine the eigenvectors which diagonalize 
the perturbation inside each degenerate subspace. We will now show that this is not the case: 
in fact, that one can work with the transformed Hamiltonian only knowing an eigenvector 
basis of the unperturbed Hamiltonian Hq, as is usually the case. 

To show this, we will explicitly diagonalize H^ in the degenerate subspace, spanned by 
\^kqa) m our example. Its matrix elements are: 

« CT I# {2) I^V> = 

= + \{(Kqa\ [Ri,h] \K q „) + <BU [Ri,h] ~ ^Im(AtJ [Ruh] \BtJ (42) 

= £t + \ [Ru h] \A+ qa ) + (B+,| [R U h] \B + kqa )} + ^Im(At qa \ [R u h] |£+ CT > (43) 



(riq«\ Hi2) \^) = 

10 



kqa I 



= \ [Ri> h] \At qa ) - (BtJ [Hi, h] \BiJ] - l -^Re(Ai qa \ [R u h\ |B+ CT > (45) 

Now, the terms in can not transform \A kqo ) into \B kqa ) because neither Ri nor 

Ii change the number of \q) phonons, this being instead the effect of [JZ kqa , {l^qa + ^kqa)]- 
Therefore, off-diagonal matrix elements (Af\ [R 1: Ii] \B k ) vanish. Conversely, the diagonal 
elements {A kq \ [R 1: h] \A kqa ), {B kq \ [Ri, h] \B kqa ) do not vanish because h] includes 
diagonal terms like tf p b p (l — n r _ pa ) n ra (1 — 5 pq 5 r k) which do not necessarily vanish in the 

\At qa )i 



\B kqa ) states. It follows that: 



kqa 



kqa I 



kqa I 



Ztqo + 7 { (4^1 ^ J l] + K*\ [*, h] \B + kqa )}) 



(46) 



kqa 



kqa I 



kqa I 



[(At q J [Ruh] \A + kqa ) - (B+J [R l7 h] \B+ qa )] 



(47) 



4 L^'-kqal I L l> 1J l^feqa/ X^fcgal L^l' l^fcg, 

Thus, diagonalization of if ( 2 )in the degenerate subspace leads to the following secular equa- 
tion: 

r Ka + (A + B)-X] (A-B) 

(A - B) [E£ qa + (A + B) - A] 



det 







where 



A=- A {At qa \[R^h]\At q(T ) 
Therefore, the eigenvalues are: 



B^\(Bt q J[Ri,h]\BtJ 



VA = Etqa + 2A = Etqa + [R U h] \A+ ga ) 

VB = Etqa + 25 = Etqa + \kKqa\ \ R ^ h] \B^) 

The eigenstate \& kqcr ) = o,\ip kqa ) + b\^ kqcr ) has coefficients set by 

[££ qa + (A + B)- VAB ]a + (A-B)b = 
11 



(48) 



(49) 



(50) 



(51) 



or, a = ±b yielding a = 1/ a/2 corresponding to the normalized eigenvectors 

i*y = ^ do + 1^>] = i4k> (*u = (avi ( 52 ) 

l*k> = ^ [l^> - l^>] = = -^ B U (53) 

Our final result is: 

<*W #e//l*W = (^l^e//l^> = ^ + ^JJ h] \At qa ) (54) 

(^i^//i^) = (^\H eff m qa ) = st q . + [^i. j i] ( 55 ) 

The conclusion is that , in evaluating the spectrum of if ( 2 We can use the eigenstates of H 
and simply ignore the "dangerous terms" [JZ kqa , (l^qcr + ^kqa)] > because only the "innocuous 
terms" will modify the non-interacting spectrum. 



V. CONCLUSIONS. 



We have explicitly shown by working out in detail an example based on the Su-Schrieffer- 
Heeger Hamiltonian, how to complete the definition of the transformation generator when 
degenerate eigenstates of the non-perturbed Hamiltonian are present. Namely, these states, 
which would apparently cause a divergence in the effective pairing between fermions and 
bosons generated by transforming "a la Frohlich" an interacting Hamiltonian, actually con- 
tribute vanishing matrix elements to the effective Hamiltonian matrix. The matrix elements 
of the latter can therefore be evaluated by using the eigenstates of the basic Hq, and ignor- 
ing the "dangerous terms". The conclusions do not depend on the size of the degenerate 
eigensubspace, as our argumentation can be straightforwardly generalized for the case of 
more than two degenerate states. Notice that, even if we have explicitated the eigenstates 
) as product of separate fermionic and bosonic eigenfunctions, actually our reasoning 
holds for any structure of \X£ ). For non-factorized eigenstates, the numerical coefficients 
in the formulae above would be different, but the conclusion would still be the same. The 
procedure has been successfully applied also to other complex fermion-boson Hamiltonians, 
e.g. the one for interacting ferromagnetic spin waves and electrons considered in 3 . Thus, 
at low temperatures renormalized excitation energies of the real interacting bosons can be 
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described by the eigenvalues of the effective second-order boson Hamiltonian, obtained by 
averaging over the fermion ground state wavefunction. 

As a final observation, notice that when we isolate the dangerous term in the perturbation 
(Eq.18) it comes weighed by a factor l/y/N. So in the infinite-lattice limit this part might be 
assumed to give in any case an irrelevant contribution . In recent times, a great deal of work 
has been done on finite or low-dimensional systems (mesoscopic, quantum dots etc.) where 
such a justification for neglecting the dangerous terms would not hold. Our demonstration, 
however, holds even for a single-particle system, because it concerns the matrix elements, 
irrespective of the kind of system where they are evaluated. 
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